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'-^ ! Electroweak baryogenesis can occur in the "adiabatic limit," in which expand- 

I I ing bubbles of true vacuum are assumed to have rather thick walls and be slowly 

qj I moving. Here the problem of calculating the baryon asymmetry in this limit is 

i-^ ' reconsidered. A simple prescription for obtaining the relevant kinetic equations 

is given. An additional suppression beyond that usually assumed is found. This 

^ ■ arises because the generation of an asymmetry requires violation of approximately 

conserved currents by Higgs expectation values, which are small near the front of 
the walls. As applications, the asymmetries in multi-Higgs models and the minimal 
supersymmetric standard model are estimated to be proportional to a^ and a^ 
respectively. Also, in this limit the baryon asymmetry in the Minimal Standard 
Model is extraordinarily small. 
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1. Introduction 

The origin of the baryon asymmetry of the universe is a long-standing prob- 
lem. With the recognition that baryon number is badly violated in the standard 
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model at high temperatures, there have been many proposals for how an asym- 
metry might be generated at the electroweak phase transition. Virtually all of 
these proposals assume that the electroweak phase transition is first order, pro- 
viding the required departure from equilibrium. CP violation already exists in 
the minimal standard model (MSM); extensions of the standard model generically 
contain further sources of CP violation. Since all three generations of quarks must 
be involved, any asymmetry which arises in the MSM should be extremely small. 
Most studies of electroweak baryogenesis have involved extensions of the MSM. 
However, there have been a number of attempts to explain the asymmetry within 
the framework of the MSM, the most recent being that of ref. 4. 

The electroweak transition, if first order, proceeds by formation of bubbles. The 
baryon asymmetry is produced in or near the bubble walls, which are the sites of 
the most significant departures from equilibrium. The various schemes considered 
to date for producing an asymmetry divide into two classes. The first is referred to 
as the "adiabatic case." If a wall is slowly moving and thick, its passage through the 
plasma is nearly adiabatic, in the sense that all quantities (but the baryon number) 
should be close to equilibrium. In this case, the time- varying Higgs field can act to 
bias the baryon-violating processes."'' ' In the second scheme, the wall is assumed 
to be thin and rapidly moving, so more significant departures from equilibrium can 
occur. Cohen, Kaplan and Nelson have pointed out that in this limit, scattering 
of particles from the wall can produce an asymmetry in front of the wall for some 
(approximately conserved) quantum number. This, in turn, can bias the baryon 
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violating processes. Other non-adiabatic scenarios have been considered in refs. 
9 and 10. 

In this paper, certain aspects of the adiabatic case will be considered. The 
adiabatic picture has been most carefully developed in refs. 6 and 7; it is closely 



related to the "spontaneous baryogenesis" scheme of Cohen and Kaplan. In ref. 
6, the problem was analyzed by obtaining the leading coupling between the time- 
varying Higgs field and Chern-Simons number. In order to make the analysis as 
simple as possible, models were considered in which this coupling arose from loops 
of heavy fields (massive compared to the temperature of the transition). However, 
while the analysis is simple in this limit, the resulting asymmetry is quite small, for 
two reasons. First, there is simply the suppression by the masses of heavy particles. 
But in addition, there is a suppression by four powers of coupling constant, g. This 
latter suppression results because the operator which gives rise to the coupling is 
quadratic in gcj) (where g denotes a gauge or Yukawa coupling and denotes the 
Higgs field). The baryon number violating process cuts off at a value of the scalar 
field, (pco (in the notation of ref. 3), for which 

g<Pco ~ (^wT, (1.1) 

i.e. for a rather small value of the Higgs field. There has been some debate in 
the literature as to how large, numerically, this suppression may be. This question 
will be taken up later; here the parametric dependence on the couplings, g^ will be 
determined. 

In many models the coupling of the scalar field to the Chern-Simons number 
arises due to light fields. While this case is potentially more promising, it is inher- 
ently more complicated, involving all the subtleties of real-time, finite temperature 
field theory. In order to deal with this, an attractive method has been proposed by 
Cohen, Kaplan and Nelson in ref. 7. The results for the asymmetry are distinctly 
more promising; in particular, it is not obvious from the analysis of ref. 7 that the 
suppression by powers of coupling constant mentioned above is obtained. 

However, this method has limitations. In the multi-Higgs model treated by 
many authors, for example, it cannot be valid for very small quark mass, since a 
non-zero asymmetry results even as the top quark mass tends to zero. This fact 
has already been remarked by the authors of ref. 7, who argue that for small 
mass there will be suppression by powers of Yukawa couplings. But a non-zero 



asymmetry is also obtained for vanishing Higgs expectation value, implying an 
additional suppression proportional to powers of ipco has been neglected. The 
method also does not lend itself to the treatment of a number of other interesting 
situations, such as baryogenesis in the MSM. 

In this note, baryogenesis in this adiabatic, light field, case will be reexamined. 
In this discussion, (p^o will be assumed small. While this leads inevitably to rather 
small asymmetries, it significantly simplifies the analysis. This is because the B- 
violating interactions switch off rapidly, in a time of order 
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where d, is the wall thickness and v its velocity, and A0 is the total change in 
across the wall (typically of order T). Even if the wall is rather thick, tco ~ T~^. 
During this time, the particle number densities can change only by small amounts. 
The kinetic equations can then be written down and solved almost trivially. 

The case where the cutoff is not so small can be analyzed by a more complete 
set of kinetic equations then will be considered below. This will be postponed for 
a future publication. 

Given the assumption of small (pco-, we first note that a more direct, naive 
treatment of the problem leads to the sort of suppression by powers of (pco observed 
in the large mass case. This naive treatment involves two steps. The coupling of the 
time-varying field to the Chern-Simons number is first computed. A rate equation 
for the baryon number is then developed by determining how this coupling biases 
the baryon-violating rate. We next consider the approach of ref. 7, in which 
the theory is rewritten so as to exhibit the coupling of the Higgs field to certain 
currents. Here the issue will be to derive suitable rate equations. In fact, the 
results of the naive treatment are recovered. 

In the end, then, the asymmetry in the adiabatic limit for the multi-Higgs 
model is not of order a^, but rather a^ (times CP- violating phases and dynamical 



factors involving sphaleron rates, the bubble profile and velocity, etc). While 
this sounds alarmingly small, it may yet correspond to an acceptable asymmetry. 
For example, general arguments give that the high temperature baryon-violating 
rate goes as naf^, but inclusion (or not) of factors of Att is currently a matter of 
prejudice. Indeed, the only calculation of this rate is that of ref. 13; the authors of 
this paper have recently expressed skepticism as to the validity of their result. As 
for the significance of the suppression, some authors have argued, by interpolating 
a formula valid for large 0, that 

g(f)co ~ 7a^ ~ 0.2 T. (1.3) 

One of us has argued that formulas which interpolate between the small and large 
limits are likely to give 2 — 3 instead of 7 in eq. (1.3). ^ The same estimate, 
however, tends to give quite a large value of k. At best, these estimates are just 
educated guesses. More extensive simulations are essential to give results reliable 
even as to order of magnitude. 

Apart from revising earlier estimates in certain models, the methods developed 
here may be applied to the MSM. Recently it has been suggested that the MSM 
may lead to a suitable baryon asymmetry. Perturbative treatments of the phase 
transition suggest that the wall is indeed thick and slowly moving, so that the 
adiabatic treatment should not be unreasonable. In this case, the asymmetry 
turns out to be extremely small (n^/s ~ 10^^^ or so). In ref. 4, however, it was 
argued that perturbation theory is not a reliable guide to the physics of the phase 
transition, and that the wall is fast and thin. We will comment on this possibility, 
but are not in a position to make any definitive statements. 



• Here a^ denotes a generic coupling; in the niulti-Higgs model, for example, this would be 
multiplied by {Xt/gwY- 



2. Naive Treatment of the Multi-Higgs Model 

For definiteness, following ref. 7, a multi-Higgs model with a CP violating 
Higgs potential will be considered. Only one Higgs, 0i, couples to quarks (to avoid 
flavor- changing neutral currents). In the bubble wall, 0i = pi{x)e^^^^' , where 9i{x) 
arises from the CP violation. The coupling of 6i to FF may be calculated from the 
triangle diagram of flg. 1. The propagators in the loop are understood as suitable 
finite-temperature Green's functions. The real-time expression can be obtained 
by analytic continuation of the Euclidean time result. In the present case, this is 
straightforward. The diagram is perfectly finite in the infrared and ultraviolet, and 
for slowly varying 9i (compared to T^^) can be expanded in a power series in the 
momenta. The result is necessarily quadratic in rrit = XtPi because of the required 
chirality fiip {rrit is the effective t quark mass. At is the t quark Yukawa coupling). 
The calculation is quite straightforward (related calculations, for example, have 
been performed in ref. 16) and one obtains 

Ce = ae^^^,-,o(^Y (2.1) 

T2 327r2 \T J ^ ^ 

where a = jj^Ci'^)- 

How might this coupling bias the sphaleron process? Integrating by parts, Cg 
may be written in terms of the Chern-Simons number 



Ce = a [j;^ \ doOmcs. (2.2) 

Suppose that the coefficient of Ucs here is very slowly changing with time. This 
gives, effectively, a chemical potential for Chern-Simons number, 

Pes = a {ml/T^) dodi. (2.3) 

In an (anti-)sphaleron transition, rics changes by (-1)1. Suppose F is the rate for 
transitions changing rics at equilibrium. Considerations of detailed balance then 



* Ref. 16 indeed calculates many of the couplings necessary for the analysis. However, they 
use the results of ref. 7 without modification to obtain rate equations, and thus their final 
results differ from ours. 



give that, when the particle densities are small compared with the equilibrium 
values (as is the case for our assumption of short times), the difference in rates for 
transitions changing rics by +1 and —1 is 

^^.c,/3r (2.4) 

The corresponding change in baryon number is three times larger. 

This result will be obvious to many readers, but a brief description of the 
derivation is perhaps useful. Near equilibrium, for small number densities, states 
differing in rics by one unit differ in free energy by an amount ^cs- Yet at equi- 
librium, the rate of transitions increasing and decreasing the free energy must be 
equal. Thus the ratio of these rates must equal the ratio of Boltzmann factors for 
the two states, 

e"^'^. (2.5) 

For small /i, this just gives eq.. (2.4). 

Since /x oc ml oc p^, this treatment should give a baryon number proportional 
to 0gQ, i.e. the small sort of rate discussed earlier. To make an estimate, we make 
the simplifying assumption of F = F(0). We then make the further simplification 
of taking F = /t(a;wT)'^ for < (pco, and F = for > (pco- This gives 

ub ~ 3aKa^\^(f)^gTA9i. (2.6) 

Here, A^i is the value of the CP- violating phase when the baryon- violating process 
turns off. The baryon to entropy ratio is of order 

It should be noted that in the multi-Higgs case, A6'i, is itself of order ((pco/T)'^ 
(times coupling constants), since in the absence of quartic couplings the Higgs 
potential is CP-conserving. 



The worst-case scenario here, in which /t ~ 1 and gipco ~ O-wT^ gives an 
unacceptably small result, of order 10~^^ for ns/ s. As remarked above, however, 
it is conceivable that /t is large, and (pco is 3 to 7 times as large, so the final result 
could well be large enough, provided CP-violating phases are large. 

Finally, we can ask how other processes affect the final asymmetry. For ex- 
ample, processes involving scattering of top quarks and Higgs fields, and QCD 
sphaleron processes, will change the numbers of left and right-handed fields. Each 
fermion species obeys a rate equation of the form 

— ^^/5/i,,rfir + 7i (2.8) 

at 

Here the first term represents the sphaleron process; (ij = for SU{2) singlets, 

di = 1/2 for doublets. As before, F is the sphaleron rate, while 7j denotes other 

processes which change the number densities. If 7j can be neglected, then summing 

over the individual rates weighted with the baryon number, eq. (2.8) reproduces 

eq. (2.4). 

The estimate described above will be valid provided that 7^ and rii are not too 
large. In order to understand this criterion, consider a particular process which 
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has been discussed recently: strong sphalerons in the two Higgs model. In this 
case, the associated terms in the rate equations are calculated just as for the weak 
sphalerons. The triangle diagram with the W bosons replaced by gluons gives rise 
to a chemical potential for "strong Chern-Simons number", n^g precisely as for rics- 
In this case, one obtains in the effective lagrangian 

CI = a^9^^^ + O (^Y (2.9) 

where a"^ = — i4C(3)- In other words, a chemical potential for the QCD Chern- 
Simons number results which is three times as large as that for SU{2). The rate 
equation for small densities is now 

-^ ~ licsm^ - 3 ha) (2.10) 

where /ij = 1/3 for color triplets and anti-triplets, hi = for leptons, and the rate. 



7 = /€^(q;sT) . In this equation, terms on the right hand side hnear in the densities 
have been neglected; this is correct provided none of the densities are of order 
HcsT"^- This equation can be integrated, as before, up to times for which = (pco- 
As argued above this corresponds to a time, tco ~ T~^ (eq. (1.2)). Apparently 
the strong sphaleron term can be neglected unless k,^ is extremely large, of order 
10^ or so (the precise value depending also on ^.) The earlier estimate of the 
baryon number is then unmodified. Other processes may be treated similarly, such 
as top quark scattering from Higgs bosons. Again, provided the estimate for tco is 
reasonable, there is no effect. The case with more complete rate equations will be 
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discussed in a subsequent publication. 

3. Spontaneous Baryogenesis 

Ref. 7 suggested a different treatment, which avoids considering directly the 
coupling to Chern-Simons number. These authors perform an anomaly-free re- 
definition of the fermion fields which eliminates the phase from the fermion mass 
terms. For example, consider again the multi-Higgs model, in the version where 
only one Higgs field couples to ordinary quarks and leptons (0i). Transforming 
each fermion by a phase proportional to its hypercharge ("fermionic hypercharge" , 
Y) eliminates the phase, 6i, from the Yukawa couplings at the price of a coupling 

0^4- (3-1) 

This appears to have induced a chemical potential for fermionic hypercharge, and 
to bias the sphaleron process. Suppose the fermionic part of the free energy is 
minimized subject to the constraints of charge conservation and separate B — L 
conservation, and including this chemical potential. The minimum lies at a nonzero 
value of the baryon number. Detailed balance arguments similar to those given 
above yield an equation for the rate of change of baryon number. Even without 
writing this equation however, it is clearly not the one above. In particular, it does 
not involve the modulus of the Higgs field, pi, at all. 



To understand this question better, focus again on short times. In this hmit, as 
before, the non-zero densities of various species may be neglected in writing kinetic 
equations. But in this hmit, it is clear that if the (small) Higgs vev is ignored, 
a chemical potential for hypercharge cannot bias the sphaleron process, since the 
sphaleron process does not violate hypercharge. Moreover, scattering of top quarks 
from the Higgs particles in the plasma (as suggested in refs. 3 and 7) cannot 
help, even if the scattering is rapid. The problem is that the field redefinition 
by fermionic hypercharge, while removing phases from the fermion mass terms, 
induces phases in the Yukawa couplings of the fermions to the fluctuating part of 
the Higgs field. To avoid these, write the Higgs field as 

01 = (pi + 01)6^^1 (3.2) 

where (p'l represents the (complex) fluctuating field. Neglecting pi, the lagrangian 
in terms of these fields contains the phase, 6i, only in the coupling 

where jy represents the full hypercharge current, including the scalar parts. But 
this current is conserved in any process (in the limit of small p); the chemical 
potential has no effect. 

In order to obtain any asymmetry at small times terms involving pi must be 
included. In this case, the results of the naive analysis are recovered. In particular, 
it is no longer true that jy is conserved; instead, at tree level 

df^jy = mtii-f^t + ... (3.3) 

In order to understand how much this is violated in the presence of background 
gauge fields consider, again, a finite-temperature Feynman diagram. The calcula- 
tion is identical to that encountered earlier, and gives 

We interpret this as meaning, on average, the violation of hypercharge in a sphaleron 
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transition is 

ami 



T2 



Aries- (3.5) 



So the change in free energy, on average, in a sphaleron transition is precisely that 
encountered in the naive treatment. The rate equation obtained is thus identical 
Although the discussion give here is for the particular hypercharge field redefinition 
of ref. 7, the results are more generally applicable to any model. 



4. Spontaneous Baryogenesis in the MSM and the MSSM 

A naive, perturbative treatment of the electroweak phase transition in the 
MSM, for Higgs masses larger 25 GeV or so (well below the current LEP limits, of 
course), gives a bubble wall which is thick and rather slowly moving, and therefore 
in the adiabatic regime. Having acquired some confidence in our understand- 

ing of spontaneous electroweak baryogenesis, the case of the minimal standard 
model may be considered. In the spirit of the naive treatment, the coupling of 
the Higgs field to Chern-Simons number must be found. Assuming g(f)co ~ cuT, 
the leading operator should be one with a minimal number of external Higgs fields 
(each additional loop costs roughly a factor of a, whereas a pair of external Higgs 
fields costs a factor of a^). The simplest such operator is 

CusM = ^FF. (4.1) 

Such a coupling is, of course, CP- violating, and so arise only at high orders. Indeed, 

f211 

in a manner similar to the calculation of 6 renormalization in the standard model, 
such a coupling cannot occur before 7-loop order. Six loops are required to obtain 
the Jarlskog invariant. A seventh loop involving a hypercharge gauge boson is also 
needed. These diagrams are perfectly finite in the infrared and ultraviolet, and for 
momenta small compared to T, represent the first term in a Taylor expansion of 
the amplitude in powers of the momentum. Thus no difficulty with the analytic 
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continuation of the result is expected, and we estimate 



^ ~ nalJ/g* (4.2) 



where 



J = Im Det f A^/AJy , Xjj 



aU ~io-^' 



and Xjj and Xjj are the up and down type quark Yukawa matrices. We have 
not attempted to include further suppression factors of tt, etc., since the answer is 
already extremely small (one might guess that these will be at least (27r)^^ ~ 10^^, 
since this is effectively a 21-dimensional Feynman integration). So even if k is quite 
large, and the suppression described earlier is small, the asymmetry is extremely 
tiny; 10""^^ is probably a quite conservative estimate. 

Shaposhnikov and collaborators have recently argued that the phase transition 
is much more strongly first order than perturbation theory suggests. This is an 
important ingredient in the analysis of ref. 4, where it is assumed that the appro- 
priate limit is the "thin-wall," highly non-adiabatic situation, in which scattering 
of particles from the wall is the most important process. We are rather skeptical of 
this claim. Moreover, a crude estimate of the mean free path for scattering of top 
quarks passing through the wall, gives a result of order T^\ so that even for an 
extremely thin wall, the scattering treatment may not be appropriate. However, 
it is a crucial assumption in the work of ref. 4 that perturbation theory is an 
extremely poor guide for all significant questions about the phase transition. This 
may be the case (though recent studies of two loop thermal effects have yielded only 
modest corrections ) . Still, it would be quite amazing if these non-perturbative 
effects change the final asymmetry by 22 orders of magnitude! (For a critique of 
the calculation of ref. 4 see ref. 24.) 

We close this section by extending the earlier estimates to the minimal super- 
symmetric standard model. The scale of superparticle masses is assumed to be of 
order T. The principle sources of CP violation are assumed to lie in phases of the 
/i term and soft breaking mass terms of squarks and gauginos. The estimate is 
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very similar to that of ref. 6. Diagrams as in fig. 2 give rise to couplings such as 

2 . . H1H2 FF 
Lcs = ag sm(5 ^^ ^^^^ <^^-3) 

where S represents some combination of CP-violating phases. Repeating the earlier 
estimates, we obtain for the asymmetry 

^^ . sin 5 fl^V^^Vx (8x10-8) (4.4) 



r-^ 



Here there is an additional difficulty. From limits on electric dipole moments (edms) 
of atoms and the neutron it is expected that sin 5 < 10^^ — lO""^. '' Again one 
must be fortunate with k, and (pco in order to obtain an acceptable asymmetry. Still, 
without further simulations, it is hard to rule out this possibility. A supersymmetry 
aficionado might even view this estimate as tantalizingly close, and suggestive that 
edms should not be far below the current bounds. '' 

5. Conclusions 

From all this, we conclude that electroweak baryogenesis, in the adiabatic limit, 
is less efficient than has been widely believed. Due to the existence of approximately 
conserved quantities, violated by Higgs expectation values, the asymmetries which 
arise in the adiabatic limit involve additional powers of couplings beyond those 
usually assumed. On the other hand, while formal arguments can be given to 
determine the parametric dependence on couplings, there are enormous uncertain- 
ties in the final numerical results. Even obtaining order of magnitude estimates 
requires knowledge of the sphaleron transition rates for small or vanishing Higgs 
field. These are, by definition, not accessible to semiclassical treatment. Crude 
estimates give widely varying answers, and existing calculations are, according to 
their own authors, unreliable even as to order of magnitude. Clearly, improved 
simulations are necessary if we are to know whether acceptable asymmetries can 
be obtained in the adiabatic limit. 
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In the non-adiabatic, thin wall, regime, there is good reason to believe that 
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acceptable asymmetries can be obtained. For these, k does not need to be large, 
and (f)co can be small, since the asymmetry is produced in front of the wall where 
the Higgs field essentially vanishes (indeed, this issue was stressed in ref. 8). The 
MSM and MSSM are probably far from this regime, but multi-Higgs models may 
well yield sufficiently violent transitions. If the baryon asymmetry were produced 
at the electroweak phase transition, and if the "worst case" scenario of small k, and 
(pco is correct, a significant step beyond currently popular theoretical ideas is likely 
required. 

We would like to thank A. Cohen and D. Kaplan for comments on an early 
version of this manuscript. 
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FIGURE CAPTIONS 

1) Leading diagram which couples the phase, 6i, to the Chern-Simons number. 

2) Leading diagram in the MSSM which couples the Higgs fields to the Chern- 
Simons number. 
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